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£NJ . Abstract 

i— i ■ We provide a new, simple general proof of the formulas giving the infinite sums <r(p, v) 

i-C \ of the inverse even powers 2p of the zeros of the regular Bessel functions J v (£), as 

functions of v. We also give and prove a general formula for certain linear combinations 
r~| | of these sums, which can be used to derive the formulas for <j(j>, v) by purely linear- 

algebraic means, in principle for arbitrarily large powers. We prove that these sums are 
always given by a ratio of two polynomials on v 1 with integer coefficients. We complete 
i n the set of known formulas for the smaller values of p, extend it to p = 9, and point out 

a connection with the Riemann zeta function, which allows us to calculate some of its 
values. 
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^ ■ 1 Introduction 

p: 

' In boundary value problems involving the diffusion equation the following infinite sums 

sometimes appear, 

m ; 

t-H I OO j 

most often for p = 1, where are the positions of the zeros located away from the origin 
of the regular cylindrical Bessel function J u (£), with real v > and integer p > 0. The 
sums are convergent for p > 1. As we will show in what follows, all these sums have the 
property that they are given by the ratio of two polynomials on v with integer coefficients. 
The simplest and most common example is 

In a few cases the exact expression of these polynomials are available in the literature [1] . 
The known cases are those obtained by Rayleigh, extending investigations by Euler, for 
p = 1 through p = 5, and one discovered by Cayley, for p = 8. The cases p = 6 and 
p = 7 seem not to be generally known, and will be given explicitly further on. The known 
cases were obtained in a case-by-case fashion, using the expression of the Bessel functions 
as infinite products involving its zeros. 

In this paper we will provide a simple, independent proof of all the known formulas, 
and will present a general formula from which the specific formulas can be derived, for any 
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given strictly positive integer value of p, by purely algebraic means. The proof will rely 
entirely on the general properties of analytical functions and on the well-known properties 
of the functions Jv{£,), which are generally available in the literature, for example in j2]. 



2 Definition of the Elements Involved 

The proof will be based on the singularity structure of the following analytical function in 
the complex-^ plane, 

where for the time being we may consider that v > and p > are real numbers. In case 
any of the functions involved have branching points at £ = 0, we consider the cuts to be 
over the negative real semi-axis. Preliminary to the proof, it will be necessary to establish 
a few properties of this function. 

We will consider the contour integral of this function over the circuit on the complex-^ 
plane shown in Fig. [H in the R — > oo limit. Since this circuit goes through the origin £ = 0, 
where the function will be seen to have a simple pole, we will adopt for the integral the 
principal value of Cauchy. The limit R — > oo will be taken in a discrete way, in order to 
avoid going through the other singularities of the function, which are located at £ = £ u k- 
We will see that, for large values of R and j, it is possible to adopt for R the values given 
by 

2u + l 

R = TT^— +J1T, (4) 

where for each k there is a value of the integer j such that R is strictly within the interval 
(£vfc> £i/(fc+i))- I n this way each step in the discrete R — > oo limit will correspond to a partial 
sum of the infinite sums involved. 

The proof of the expressions for cr(p, v) consists of two parts: first, the proof that the 
integral of f(p, u, £) over the circuit is zero in the R — > oo limit, for all p > and all v > 0; 
second, the use of the residue theorem. This will result in a general formula from which the 
expressions for the sums a(p, u) can be derived. We will also present a partial solution of 
the problem of deriving the formulas for a(p, v), which will take the form of another general 
formula from which these expressions can be derived algebraically. 

3 Properties of f(p,v,£) 

Let us establish a few important properties of f{p, u, £), starting by its behavior under the 
inversion of the sign of £. We start with the analogous property of J v (£), for which we have, 
using the Maclaurin series for these functions [3], which converges over the whole complex 
plane, 

M-0 = (-irMO- (5) 

Using this in the expression for f(p, v, £) we get 

/(p,f,-0 = -/(P,^0. (6) 
that is, f(p, v, £) is an odd function of £, for all p and all v. 
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Figure 1: The integration contour in the complex-^ plane, showing the various parts of the 
circuit C and some of the singularities of f(p,v,£)- 



Next we show that f(p,v,£) has a simple pole at £ = 0. Since J u+P (^), and 
J u {i) are analytical functions over the whole complex-^ plane, it follows that f{p,v,C) is 
analytical over the whole plane except for those points where the denominator vanishes, 
where it has poles. These are the origin £ = and the zeros of the Bessel function 
in the denominator. Note that while for non- integer v and p the functions involved have 
branching points at £ = 0, the function f(p, u, £) never does. In order to determine the 
residue of f(p, u, £) at £ = we consider the limit 

r (p,u) = lim£/(p,i/,£) 

I>+D 



2pr(i/+p + i)' 

where T(z) is the gamma function and we used once more the Maclaurin series for J u (£). 
Since the limit is finite and non-zero, it follows that f(p, p, £) has a simple pole at £ = 0, 
and that ro(p, v) is the corresponding residue. Turning to the poles at since is 
a simple zero of J u (0i a ^ which its derivative J^(£) is different from zero, it follows that 
f(p,v,£) has a Taylor expansion around this point, with the form 

oo 

mo = (£ - iu k )Mu) +^2^ - ( 8 ) 

for certain finite coefficients c%. In order to determine the residue of f(p, v, £) at £ u k we 
consider then the limit 
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rk(p, v) 



lim (f - Zvk)f(p,v,€) 



0) 



where we used this Taylor expansion. Since the derivative is finite and non-zero at this 
limit also is finite and non-zero, and hence it follows that f(p,v,£) has a simple pole at £ u k, 
and that rk(p,u) is the corresponding residue. We can simplify this expression using the 
well-known identity [1] 



which applied at £ = since J v (£vk) = 0, results in 



(10) 



J'viZvk) = -Ju+l(Cuk)- 

It follows therefore that we have for the residues of the poles at £ u k, 



nip, v) 



Jv+p{£,vk ) 



£lk lj v+^vk) 



(12) 



We will now establish the behavior of the absolute value of the ratio of Bessel functions 
which is contained in the expression of f(p, v, £) in Eq. ([3]), for large values of R. In order 
to do this we use the asymptotic expansion of the Bessel functions [5], valid in the whole 
complex plane so long as arg(£) ^ ±7r, written in terms of R = |£| and 8 = arg(£), to the 
lowest orders, and with the trigonometric functions expressed as complex exponentials, 



MO 



2ttR 



1 + 



R 2 e* 29 



D -Rsin{e) ia{Rfi,v) 



+ e 



Rsm(6) -ia(R,9,u) 



+ 



+1 



Av 2 -l , K a {y,Z) 

R 3 e tW 



8Re 



1 6 



+ 



e -Rsin(0) e ia(R,9,v) _ e -Rsin(6») -ia(R,6,v) 



(13) 



where lZ c (v, £) and 1Z s {v, £) are certain limited functions of £ and a(R, 8, v) is a certain real 
number, given by 



a(R,0,v) = Rcos(8)-tt 



2v + 1 



(14) 



The behavior of the expression in Eq. ([13)1 for large values of R depends on the sign of 8, 
and the particular case 8 = has to be examined separately. In this particular case we 
have 



1 + 



R 2 



cos[a(R, 0, z/)] + 



4i/ 2 -l ft a (i/,0 



8i? 



+ 



ft 3 



sin[a(i?, 0, z/)] > , 



(15) 



where all the functions involved are now limited, so that for large values of R we have for 
the dominant part of </„(£), 
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ttR 



cos[a(i?, 0, v)\. 



(16) 



Note now that the points where cos[a(R, 0, u)} = are the zeros of Jj,(£), expressed in the 
asymptotic limit. We will now choose a way to take the R — > oo limit such that these zeros 
are avoided. We may simply chose for the passage of the circuit across the real axis that 
point between two zeros where cos[a(R, 0, v)] = ±1 and sin[a(R, 0, v)] = 0. Since for 8 = 
we have 



a{R,0,is) = R-ir 



2u + l 



(17) 



and we must have a(R, 0, u) = jn for some integer j, we conclude that Eq. holds, which 
will cause the crossing of the circuit and the real axis to avoid the zeros. This defines 
the R —¥ oo limit in full detail. It follows that for our purposes here we may write the 
asymptotic expansion in the case 9 = as 



MO 



ttR 



1 + 



R 2 



In the case > we put the dominant real exponential in evidence and obtain 



MO 



£ e -Rsin(0) 

2irR 



1 + 



t28 



R 2 e 



-2Rsm(9) ia(Rfi,v) 



+ e 



+ 



8Re* 9 R 3 e* 3e 



(18) 



e -2i?sin(6>) ia(R,9,v) _ -%a{R,0,v) 



(19) 



where all the functions within the brackets are now limited or go to zero in the R — > oo 
limit. Finally, we do the same thing for the case 9 < 0, obtaining 



MO 



2ttR 



-Hsin(e) 



+1 



1 + 



KMO 

R 2 e z29 



s ia(R,6,u) i 2_Rsin(6») -ia(R,9,v) 



+ 



4^-1 n 8 (v,Q 

8Re ld R 3 e l3d 



e %a{Rfi,v) _ e 2Rsin{8) -ia(R,9,v) 



(20) 



where once more all the functions within the brackets are now limited or go to zero in the 
R — > oo limit. We are now in a position to analyze the behavior of the absolute value of 
the ratio of two Bessel functions which appears in the definition of f{p,v,0- The factors 
which do not depend on v are common to the numerator and denominator, and cancel out. 
In the case 8 = we get 



Jv+piO 



MO 



i + 



R 2 



R 2 



(21) 
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so that in the R — s- oo limit we get 



lim 

R— ►oo 



MO 



1. 



(22) 



It is not difficult to verify that for both the case 9 > and the case 6* < we get this same 
value for this limit. We see therefore that the R — > oo limit of the absolute value of this 
ratio is simply 1, for all values of 6 in (— n, n). 



4 Evaluation of the Integral 

Let us consider now the proof that the integral is zero. In order to do this we will separate 
the circuit in sections and prove the result for each section. The complete circuit C consists 
of two straight sections C\ e C3, of the great semicircle C2 and of two small semicircles Cl 
e Cr of radius e around the point £ = 0. 

For the pair of straight lines C\ e C3, where we have d£ = idy with £ = x + ty, taking 
into account the orientation, we may write 



Id+Cs = / f{p,v,t)d£ 
JC1+C3 

= * f f(p,v,Z)dy + i I f(p,v,0dy. (23) 

JR J-e 



Making in the second integral the transformation of variables £ — > — £, which implies x — > 
—x e y — > —y, and since f(p, v,£) is odd, we have 

kh+c z = 0. (24) 

We see therefore that this part of the integral vanishes exactly, independently of the values 
of R and e. We are therefore free to take limits involving R or e during the calculation of 
the other sections of the integral, without affecting this result. 

Next we consider the two semicircles of radius e. We will denote this part of the integral, 
to be calculated according to the criterion of the principal value of Cauchy, as 

Ic L+ c R = [ /(P,f,0de- (25) 
Jc L +c R 

Since the function f(p, z/,£) has a simple pole at £ = 0, and is also odd, it can be expressed 
as a Laurent series around this point, with the form 

/(P,^) = + f>^+\ (26) 

^ j=0 

where ro(p, v) is the residue of the function at this point, and Cj are certain finite coefficients. 
The series is convergent so long as e is smaller than the first zero ^ v \. The sum of positive 
powers represents an analytical function around £ = 0, and is therefore regular within the 
circle of radius e. It follows that the integral of this regular part goes to zero in the limit 
£—>•(), since in this limit both each individual term of the sum and the measure of the 
domain of integration vanish. 

It follows that only the integral of the term containing the pole can remain different 
from zero in the e —¥ limit. We will therefore calculate this integral in polar coordinates, 
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with £ = eexp(i9) and (i£ = ieexp(i9)d9. Since we must use here the Cauchy principal 
value we have for this part AIc L +c R of the integral In L+ c R , 

1 f r (p,v) | 1 I r {p,v) 

2 In, £ 2 



AL 



Cl+C r 



>C L C 
2 7^/2 

0. 



ir (p, v) 



-it/2 



d9 



tt/2 



(27) 



Therefore, this part of the integral Ic L +c R also vanishes, and hence the integral Ic L +c R 
vanishes in the limit e — > 0. Since during this deformation of the circuit no singularities of 
the function are crossed, and hence the integral does not change, if follows that the integral 
is zero for all values of e smaller than £„i. 

The last section of the circuit we must consider is Ci- In this case the integral is not 
zero for finite values of R, but we may show that it goes to zero in the limit R — > oo, 
subject to the condition that for large values of R we have that Eq. ([!} holds, so that the 
circuit does not go over any of the singularities at the points £ v &. Using once more polar 
coordinates, in this section of the circuit we have d£ = i Rexp(i0)d0, where £ = Rexp(iO), 
so that the integral is given by 



Ic 2 



c 2 



n/2 



(28) 



(29) 



= iR I d6e tti f(p,v,0- 

J -tt/2 

Taking the absolute value of the integral and using the triangle inequalities we have 

I <R f 2 d6\f(p,v,0\, 

J -tt/2 

for any value of R, and hence also in the R — > oo limit. We must now consider the behavior 
of the absolute value of f(p, v, £) for large values of R. In order to do this we calculate the 
limit 

R l" W/2 d9\f(p,v,0\ 

I-tt/2 
r-TT/2 

R I dO 

-tt/2 



ft | 



< lim 

R— >oo 



lim 

R^oo 



1 

Rpl 



(30) 



As we established before, the limit of the absolute value of the ratio of the two Bessel 
functions is 1. As a consequence of this, we have for the integral over the section C2 of the 
circuit, in the R — > 00 limit, 

(•tt/2 

lim 



lift 



< 



R^oo 

0. 



tt/2 



RP 



since we have p > 0. This implies, of course, that Ic 2 = in the R 



(31) 

00 limit. We see 



therefore that the integral of f(p, u, £) over the circuit C, in the R — > 00 limit, vanishes in 
all sections of the circuit, and hence that the integral is zero in the R —¥ 00 limit, 



lim i> /(p, I /,£)d£ = 0. 



(32) 



c 
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5 Using the Residue Theorem 



Considering that in the R — > oo limit the poles with residues v) are all that exist 
strictly within the circuit, that the pole with residue ro(p,v) is the only one located over 
the circuit, and that the integral is defined as the Cauchy principal value at this pole, we 
can use the residue theorem to write for the integral 



lim j) f(p, u, i) di = 2m - r (p, v)+Vr fe (p 



(33) 



.. . ") 

k=l 

On the other hand, as we saw above the integral vanishes in the R — > oo limit, and hence 
we have 

j oo 

^r {p,jy)+J2 r k{P,^) = 0. (34) 

k=l 

We have therefore the following general result involving all these residues, substituting the 
values we calculated before for each one of them, 

r(u + i) ^ i Ju+ P (U) 



This is valid for any real value of v > and for any real value of p > 0. 

6 Proof of Some Known Formulas 

Up to this point p could be any strictly positive real number. Prom now on, however, 
we have to assume that p is a strictly positive integer. In order to further simplify the 
expression obtained above, in general it will be necessary to write Ju+p{£,uk) in terms of 
Ju+i{£,uk), which can be done using the recurrence formula of the Bessel functions [6], so 
long as p is an integer. Let us examine a few of the initial cases. For p = 1 we have simply 

iy + i) i 

^TY)~^ (36) 

so that the formula for the sum that corresponds to this case is 



fc=l ^"k 



1 (37) 



2 2 (i/ + l 

were we used the properties of the gamma function, thus obtaining a polynomial on v in 
the denominator. In this way we obtain the first of the known results, and this formula is 
therefore proven, being valid for any non-negative real value of v. For p = 2 we have 

r(i/ + i) ^ i j„ +2 (e 



— -, — 7T-- ( 38 ) 



In order to simplify this expression we write the recurrence formula as 



S 



^+2(0 = ^1^^+1(0-^(0, 



(39) 



where we exchanged v for v + 1. Applying this for £ = and using once more the fact 
that J„(6,fc) = 0, we get 



2(^ + 1) 

■Ju+2{t,uk) — J ■Jv+l{C,uk), 

C,uk 



so that in this case we have 



r(f + i) ft , i 

from which it follows that the formula for the sum that corresponds to this case is 



(40) 



(41) 



a(2,u) 



CO ^ 

k=i ^ uk 



24^ + 1)2(^ + 2)- 



(42) 



We thus obtain the second known result, which is now proven, and which also has a poly- 
nomial on v in the denominator. The proof of the first two formulas is therefore quite 
straightforward. In the p = 3 case, however, something slightly different happens. In this 
case we have 



1 Jv+3(€vk) 



rfr + i) 

2 4 i>+4) ^e uk Ju+i(uv 



(43) 



and hence we must write one more version of the recurrence formula. Exchanging v for 
v + 2 in the original formula, and applying at £ = ^ , we obtain 



Jv+2,(£,vk) — -^7 " Ju+2(£,vk) - Jv+l{£vk)- 

Kuk 



(44) 



Substituting in this the solution found for the previous case, which gives us J v +2(£,vk) in 
terms of J u +i(,^k), we get 



Jv+z{t,vk) 

In this way we get in this case the result 
I> + 1) 



2 2 {v + l){v + 2) 



vk 



Jv+l{£,vk)- 



(45) 



2 4 r(z/ + 4) 



2 2 {v + l){u + 2) 



X t4 

k=l 

2 2 (v + l)(v + 2)a{3,v) -a{2,v). 



( 2 

^vk 



(46) 



We see that in this case a linear combination of the sums of two different powers of the 
zeros appears. Using the properties of the gamma function and substituting the value 
obtained previously for a(2, v), we get for <r(3, v) 



(7(3,1/) = 



2{y 



2 6 (^ + l) 3 (z/ + 2) 2 (z. + 3)' 



(47) 
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The formula for the sum that corresponds to this case is therefore 



<r(3,f) = ^7(T 
k=i ^ vk 

25( I / + l)3( I / + 2)(i/ + 3)' (48) 

where we once more have a polynomial on in the denominator. We thus obtain the third 
of the known results, and the formula for <r(3, v) is proven. 

It is clear that we can proceed in this way indefinitely, thus obtaining the formulas for 
successive values of p. In each case it is necessary to first use the recurrence formula in order 
to write J v + p (£) in terms of J v +i(£). In general the result will be a linear combination of 
sums of several distinct inverse powers of At this point the use of the general formula 
in Eq. (135H will produce an expression for the linear combination of the corresponding sums 
a(p, u). Finally, it is necessary to solve the resulting expression for the sum with the largest 
value of p so far, using for this end the results obtained previously for the other sums. In 
this way all the formulas for the sums a(p, v) can be derived successively by purely algebraic 
means, resulting every time in the ratio of two polynomials, with the one in denominator 
completely factored. 

7 Proof of a General Formula 

It is possible to systematize the resolution process described above to the point where a 
general formula for the linear combination of the sums a(p, v) can be written. This is based 
on a systematization of the general formula for the ratio of Bessel functions, which is found 
to be 

■WU) = vv, )q [(p-i)-g}ir(u + p- q ) ^- 2 * 

Ju+l(U) ^ ' [(p-l)-2q}lq\T(u + q+l)\Uj ' 1 ' 

where qM = (p — l)/2 for odd p and qu = (p — 2)/2 for even p, and for which we will 
provide proof in what follows. The use of the general formula in Eq. (|35p then produces a 
corresponding general formula for the linear combination of the sums a(p, u), 

2fT(u+p+l) ^ [(p - 1) - 2q\l q\ Y[v + q + lj 

Note that the left-hand side of this equation can be written as the inverse of a polynomial 
on u, with integer coefficients, with the simple use of the properties of the gamma function. 
On the other hand, the coefficients on the right-hand side can all be written as polynomials 
on is, with integer coefficients, since we have for the arguments of the two gamma functions, 
in the numerator and in the denominator, 

(v+p-q) = (v + q + l)+p-2q-l, (51) 

where p — 2q — 1 is an integer whose minimum value is for odd p, and 1 for even p. It 
follows that, once the equation is solved for a(p,u), resulting in 
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T(v + p) _ I> + 1) 

T(is + l) a(P,V) ~ 2 2 Pl>+p+l) + 



i )g [( p _i)_ g ]ir(z. + p-^ 



9= 



22g [(p_ i) _ 2 g ]! g ! r(j/ + g + 1 



and assuming that the previous sums all have this same property, the expression for this 
sum will have the form of the ratio of two polynomials on v, with integer coefficients. Hence, 
since we saw that this is valid for the first three sums, by finite induction it is valid for all 
the sums. 

This set of equations, taken for all strictly positive integer values of p, forms an infinite 
linear system of equations in triangular form, that can be solved iteratively in order to 
obtain closed forms for a(p, v) in a purely algebraic way, in principle for arbitrary integer 
values of p, although for large values of p the algebraic work involved can be very large. 
However, it is straight, direct algebraic work, well suited for a computer-algebra approach. 

We will now prove these two general formulas. Since the general formula in Eq. (|50p 
follows from the general formula in Eq. (09]), it suffices to prove the latter. We can do this 
by finite induction. Since the upper limits of the summations involved depend on the parity 
of p, it is necessary to consider the two cases separately. The first step is to verify that 
our general formula reproduces the correct results for the first three cases, which we have 
already derived individually. Applying the general formula in Eq. (|49p for p = 1, in which 
case we have qu = 0, we obtain at once 

1 — TT^ = lj ( 53 ) 

which is obviously the correct result. Applying now the same general formula for p = 2, for 
which we also have qu = 0, we get 

1 — TT^ = ~? ' (54 ) 

which is also the correct result. Finally, applying the general formula for p = 3, in which 
case we have qu = 1, we obtain 

Ju +3 (U) _ 2> + 2)(^ + l) 
Jv+\\Kvk) k vk 

which once more is the correct result. It suffices now to use the recurrence formula of the 
Bessel functions to show that the formula for p follows from the previous formulas, for p — 1 
e p — 2. We start with the case in which p is even, and writing explicitly the upper limits 
of the sums, we have 

■Wi(U) (P y 2 r 1v , [(p-2)-g]!rfr + p-i-g) / 2_Y~ 2) ~ 2q f ~n 

Ju+l(U) ^ 1 ' [(p-2)-2q]\q\T(u + q+l) VW ' ^ ' 

J» +P -2(U) (P y 2 f n , [(p-3)-g]!r(i/ + P -2-g) / 2 \fr-3)-2g 

Ju+i(U) ^ 1 ' [(p-3)-2q}\q\r(u + q + l) \Uj ' 1 ' 

Writing now the recurrence formula which gives the function J u +p(£,uk) i n terms of Jv+p—i(£uk) 
and J v+p -2(Cuk), and substituting Eq. (|56l) and (j57]) . we get, after some manipulation of 
the indices of the sums, 
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Jv+p(£vk) — — y~ Jv+p-l(£vk) — Jv+v-2(tivk) 

Ju+p{iuk) _ 2(u + p — 1) J v+p -i(tl vk ) Jp+p-2{iuk) 



Jv+l(£vk) ivk Jv+l(£vk) Jv+l(€uk) 



[(p-l)-2q}\q\T(u + q+l) 
(v + p-x)[(p-l)-2q]+q(y + q) ( 2 \ 



[{p-l)-q)[u+(p-l)-q) VW 
It is easy to verify that we have for the second fraction in this sum 

(u+p-l)[(p-l)-2q]+q(u + q) 



(p-l)-2g 

) ■ (58) 



[(p-l)-q}[u + {p-l)-q] 

We therefore conclude that 



1. (59) 



■WU) = iP ^ /2 ( _, )q l(p-i)- q ]ir(v + P- q ) f±\ {p - 1] - 2q (m) 

Ju+l(U)~ ^ l 1 [{p-l)-2q]\q\T{u + q+l) \u) ' 1 ' 

thus proving the general formula for even p. For odd p, once more writing explicitly the 
upper limit of the sums, we start from 

Ju+p-l(U) (P y 2 r , )Q [(p-2)-q]lT(u+p-l-q) f J_\ <*-*>~* 

Jv+l&k) ^ 1 ] [{p-2)-2q]\q\T{u + q+l) VW ' 1 ' 

J» +P -2(U) (P y 2 r lV , [(P-3)-g]!r(i/ + p-2-g) / 2 \ ^ 
Ju+l(U) ^ 1 ] [(p-3)-2q]lq\T(u + q+l) \u) ' k ' 

Writing once again the recurrence formula which gives J v +p(Cuk) hi terms of J u +p~i(£,vk) an d 
Ju+p-2(S,uk), and substituting Eq. (JoTj) and (j62j) . we get, after some similar manipulation 
of the indices of the sums, 

Jv+p{ivk) = — T~~~ Jv+p—xi^vk) ~ Jv+p—2i£,vk) =^ 

C,uk 

Ju+p{iuk) _ 2{y +p - 1) J u +p-l(£,vk) Jv+p-2{ivk) 



Jv+x{fivk) ivk Ju+x(Cvk) Ji>+x(Cvk) 

(p 3^/2 [(p-l)-g]!r(t/ + p-g) 

^ 1 ' [(p-l)-2q]\q\T{u + q+l) 

{u+p-l)[(p-l)-2q]+q(u + q) ( 2 \ 



(p-X)-2g 



[(p-l)-q}[u+(p-l)-q\ VW 
+(_l)(p-i)/2. (63) 

The second fraction within the last summation, which does not involve factorials, is the 
same as before, and therefore is equal to 1. It follows that we have 

J» +P (U) iP s^\ 1V , [(p-l)-q]\T(v+p-q) / 2 



Jv+i{iuk) (/ 



V r iv [(p-i)-g]'r(^+p-g) py 

^ J [(p-1) -2?]! ?! 1^ + 9 + 1) UJ 
+(_ 1 )(p-l)/a. (64) 
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It is not difficult to verify that the additional term that we have here is in fact equal to the 
argument of the summation in the case q = (p — l)/2, so that we may merge it with the 
summation and this obtain 



Ju+i{U) ^ V ' [(jp-\)-2q]\q\T{v + q + \) \U 

thus proving the general formula in this case. This completes the proof of the general 
formula in Eq. (|49p . from which follows the general formula in Eq. (|50p for the linear 
combination of the sums a(p, u), which is therefore proven as well. 

8 Some Particular Cases 

We may now use the general formula in Eq. (|50p for a(p, u) in order to write explicitly a few 
cases which are not found in the current literature. With a little help from the free-software 
algebraic manipulation program maxima, we get the following two results, thus completing 
the sequence of known results up to p = 8, 

( v 21^ 3 + 181^ + 513^ + 473 

a{,U) 2 11 (zy + l)6(z/ + 2)3(z/ + 3) 2 ( l / + 4)(z. + 5)(^ + 6)' [ ' 

33z^ 3 + 329^ 2 + 1081t/ + 1145 

a " ,V ' ~ 2 l2 {u + \y{v + 2) 3 (i/ + 3) 2 (j/ + 4)(i/ + 5)(i/ + 6)(i/ + 7) ' ( } 

We now point out that our results for a(p, v) are valid for all real values of v, not just for 
the integers. Therefore, exchanging v for v + 1/2 we may obtain formulas that are valid 
for the zeros of the regular spherical Bessel functions j v {£), since we have the well-known 
relation 

3u(0 = <Mju+i/2(t) (68) 

between these two families of functions. In particular, using the value v = 1/2 we obtain the 
results for jo(£), whose zeros are given by kw, since this particular function is proportional 
to sin(£) [7]. In this way we obtain a direct relation between our results and the Riemann 
zeta function, for certain real integer arguments of CO 2 )- in fact, we have 



k=l 

C(2p) = tt 2 V(p,1/2). (69) 



Using the formulas we obtained here for a(p, u) in the case v = 1/2, we obtain for example 
the values 

6917T 12 

C(14) = - R =-=^ . (71) 

sv ' 3 6 x 5 2 x 7 x 11 x 13 y ' 

Finally, using our general formula for the case p = 9 and solving for cr(9, is), once more with 
some help from the free-software program maxima, we obtain 
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<t(9,i/) 



(72) 



Q 9 (uy 



where the two polynomials are given by 



715 v % + 16567 i/ 5 + 158568 v 4 + 798074 v 6 



+ 



Qg(v) 



+2217079 v 2 + 3212847 v + 1893046, 
2 17 (^ + 1)> + 2)> + 3)=V + 4) 2 x 
x(u + 5)(i/ + 6)(i/ + 7)(i/ + 8)(i/ + 9). 



(73) 



(74) 



Up to the case p = 8 it is possible and in fact fairly easy to verify the resulting formulas 
numerically with good precision, with the use of standard computational facilities. However, 
in this p = 9 case it is just too difficult to verify this formula by numerical means, except 
for v = 0, using the usual double-precision floating-point arithmetic. In order to do this one 
would have to use quadruple precision or better numerical arithmetic. The difficulty seems 
to lie in the direct numerical calculation of the sum a(9,v), not in the evaluation of the 
ratio of polynomials. Hence, the results discussed here acquire an algorithmic, numerical 
significance, enabling one to easily calculate the values of the sums. 
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10 Historical Note 

From the information supplied by Prof. Muldoon it seems that, since the time of Euler, 
aspects of the subject of these sums has been discovered and rediscovered, possibly more 
than just once. The explicit expressions of the sums as functions of v were given, for the 
first 12 cases, by Lehmer [8]. The relation of the sums to the diffusion equation was first 
established by Kapitsa [U[T0], who also worked on the calculation of the sums. The recursive 
general solution presented here seems to be equivalent to a known recursion formula, first 
derived by Meiman [11] and later rediscovered by Kishore |12| . 
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